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ABSTRACT 

We present a study of 33 Kepler planet-candidate host stars for which asteroseismic 
observations have sufhciently high signal-to-noise ratio to allow extraction of indi¬ 
vidual pulsation frequencies. We implement a new Bayesian scheme that is flexible 
in its input to process individual oscillation frequencies, combinations of them, and 
average asteroseismic parameters, and derive robust fundamental properties for these 
targets. Applying this scheme to grids of evolutionary models yields stellar properties 
with median statistical uncertainties of 1.2% (radius), 1.7% (density), 3.3% (mass), 
4.4% (distance), and 14% (age), making this the exoplanet host-star sample with the 
most precise and uniformly determined fundamental parameters to date. We assess 
the systematics from changes in the solar abundances and mixing-length parameter, 
showing that they are smaller than the statistical errors. We also determine the stel¬ 
lar properties with three other fitting algorithms and explore the systematics arising 
from using different evolution and pulsation codes, resulting in 1% in density and 
radius, and 2% and 7% in mass and age, respectively. We confirm previous findings 
of the initial helium abundance being a source of systematics comparable to our sta¬ 
tistical uncertainties, and discuss future prospects for constraining this parameter by 
combining asteroseismology and data from space missions. Finally we compare our 
derived properties with those obtained using the global average asteroseismic observ¬ 
ables along with effective temperature and metallicity, finding an excellent level of 
agreement. Owing to selection effects, our results show that the majority of the high 
signal-to-noise ratio asteroseismic Kepler host stars are older than the Sun. 

Key words: Asteroseismology — stars: evolution — stars: oscillations — stars: plan¬ 
etary systems — stars: fundamental parameters—planets and satellites: fundamental 
parameters 
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1 INTRODUCTION 

The CoRoT and Kepler missions have revolutionised the 
quality and quantity of data that are available for the analy¬ 


in the ecliptic plane have already started with K2 (Howell 


sis of solar-like oscillations (see e.g., Chaplin & Miglio 20131. 


As a result, asteroseismology of solar-type stars has ma¬ 
tured into a powerful tool to help characterise extrasolar 
planetary systems. Even a relatively straightforward analy¬ 
sis of solar-like oscillations can usually provide the surface 
gravity and mean density of an exoplanet host star with 
a typical precision much higher than the constraints avail¬ 
able from transit modelling (e.g., Huber et al. 2013b| Van 
Eylen et al. 2015, submitted). With an independent effec¬ 
tive temperature constraint, a precise stellar radius from as¬ 
teroseismology can be obtained and used to determine the 
absolute planetary radius from transit photometry. The as- 
teroseismic mass provides the absolute scale of the orbit, 
and when combined with transit-timing or radial velocity 


measurements yields the absolute planetary mass (Carter 


et al. 2012 Marcy et al. 20141. The asteroseismic age can 


be used to assess the dynamical stability of the system, and 
to establish its chronology with respect to other planetary 
systems. For the best and brightest targets, asteroseismol¬ 
ogy can also provide an independent determination of the 


stellar rotation rate and inclination (Chaplin et al. 2013 
Van Eylen et al.j2014| Davies et al.|2015l. This can be used 


to probe spin-orbit alignment for systems without Rossiter- 
McLaughlin (R-M) measurements during transit, and it pro¬ 
vides the full three-dimensional orientation when R-M mea¬ 
surements are available (e.g., [Benomar et al. 1 |20T4l|L^ 
et al. p0l4| . As a further matter, lower limits on the surface 


gravities of exoplanet-host stars can still be placed in the 
event of a non-detection of solar-like oscillations ( |Campante| 
et al.|20T4 1. 


Prior to the CoRoT and Kepler missions, solar-like os¬ 
cillations were extremely difficult to observe for anything 
but the Sun, and data only existed for a handful of the 


brightest stars in 

the sky (e.g., Kjeldsen et al.|1995| Bouchy 

& Carrier||2001| 

Carrier, Eggenberger & Bouchy||2005| |Ar- 

entoft et al.|2008 

1. With the advent of space-born asteroseis- 


mic observations, thousands of detections have been possi¬ 


ble for stars in different evolutionary stages (e.g., De Rid- 
der et al. 2009 Huber et al.||2011 1. The focus of recent ef¬ 


forts has been to develop uniform data analysis and mod¬ 


elling strategies (Appourchaux et al. 2012 Chaplin et al. 


2014), and to optimise the precision of asteroseismic infer¬ 


ences by matching the individual oscillation frequencies or 


ratios of characteristic frequency separations (Silva Aguirre 
et al.|2013 Metcalfe et al.|2014 l. These techniques typically 


improve the precision by a factor of two or more over meth¬ 
ods that only use the global oscillation parameters (e.g., 
Mathur et al.| [^012| [Silva Aguirre et al.||2013l [Lebreton fcj 
Goupil 2014), but the absolute accuracy is difficult to as¬ 


sess. Initial comparisons for the small samples of dwarfs and 
subgiants with well-determined distances and interferomet¬ 
ric observations suggest consistency at the level of 4-5% in 


stellar radii (Silva Aguirre et al. (2012); Huber et al. (2012), 


see also Miglio et al. (2013aI for a summary). 


Future space-based photometry and astrometry, as well 
as complementary data from ground-based networks, will 
soon yield improved asteroseismic constraints for many ex¬ 
oplanet host stars. Observations of bright stars and clusters 


et al. 2014), and the TESS mission will soon obtain sim¬ 


ilar data over the entire sky (Ricker et al. 2014). At the 
same time, the Gaia mission will determine improved paral¬ 
laxes for all of these targets ( Perryman et al.|2001 1, yielding 
luminosity constraints that could help break intrinsic pa¬ 
rameter correlations in stellar models. The PLATO mission 
will build upon these successes with extended time-series of 
many helds around the Galactic plane ( Rauer et al.||20l4 |, 
while ground-based spectroscopy from the Stellar Observa¬ 


tions Network Group (SONG, Grundahl et al. 2008) will 
complement the satellite data for the brightest stars in the 
sky. Now is the time to address the underlying sources of 
systematic uncertainty in stellar properties determination, 
before this next deluge of observations. Exploration of biases 
in asteroseismic models and fitting methods for individual 
targets has been made (e.g., Miglio fc Montalban|20d5 Le- 
breton fc Goupil|2014 1, and we intend in this study to com¬ 
plement these efforts using a larger sample of stars together 
with a range of evolutionary and pulsation codes, and fitting 
algorithms. 

In an accompanying paper (Davies et al. 2015, in prepa¬ 
ration) the extraction of individual frequencies, combina¬ 
tions, and derivation of correlations has been performed for 
a sample of Kepler Objects of Interest (KOIs) with high- 
quality asteroseismic observations. We make use of these 
data to determine a robust set of fundamental properties, 
including ages, by introducing a new BAyesian STellar Al¬ 
gorithm (BASTA) and grids of evolutionary tracks specially 
constructed for this purpose. Our modelling strategy ini¬ 
tially adopts this fitting algorithm and applies it to a fixed 
set of observables and modelling ingredients, and uses the 
results as the comparison basis for quantifying the impact 
of additional input physics (e.g. diffusion, convective over¬ 
shoot, mixing-length, composition) in a systematic manner. 
We further explore the effects introduced by choosing dif¬ 
ferent sets of asteroseismic observables, evolutionary and 
pulsation codes, as well as treatment of the pulsation data 
by determining stellar properties with three asteroseismic 
pipelines widely used in the literature. 

This paper is organised as follows. In Section we de¬ 
scribe the target sample and define the observational quanti¬ 
ties to be matched, while our new model grids and Bayesian 
analysis scheme are outlined in Section The main results 
are discussed in Section including a study of systematic 
uncertainties introduced by changes in the input physics 
and a thorough comparison with the results of asteroseis¬ 
mic htting algorithms described in Appendixj^ We use our 
derived stellar properties to check for consistency with the 
results from empirical scaling relations in Section and ex¬ 
plore possible correlations of planetary properties with age 
in Section Conclusions and closing remarks are given in 
Section 0 


2 TARGET SELECTION AND DATA 

Our sample has been extracted from the 77 exoplanet host 


stars presented in Huber et al. (2013b I. These stars show 


clear signatures of stochastically excited oscillations with 
power spectra characterised by a gaussian-shaped envelope 
centred on the frequency of maximum power t'max. Pulsation 
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Figure 1. Echelle diagram of two selected targets from the sample. Symbols depict the different angular-degree values £ = 0 (circles), 
£ = 1 (triangles), and £ = 2 (squares). Uncertainties in some frequencies are smaller than symbol sizes. Left: Kepler-409 shows the 
quasi-vertical alignment of its ridges as expected from pure p-mode frequencies. Right: Kepler-36 presents signature of modes with mixed 
character in the £ = 2 ridge at ~1070 /iHz. 


modes of the same angular degree £ and consecutive radial 
order n are separated in frequency by a roughly constant 
spacing, known as the large frequency separation, defined 
as 




( 1 ) 


Huber et al. 12013b I determined (for all 77 targets) the 


~0.3 /rHz (and below 0.1 /rHz in the best cases), while 
fractional uncertainties in the same range for the frequency 
ratios are ~7% (and smaller than 3% in the highest SNR 
stars, further details can be found in the accompanying 
paper of Davies et al. 2015, in preparation). The advan¬ 
tage of using the ratios instead of the individual modes 
of oscillation is that they effectively suppress the contribu- 


effective temperature and surface composition [Fe/H] tions from the poorly modelled surface layers (see e.g., Rox- 


from high-resolution spectroscopy, as well as the global av- 

burgh|2005||0tf I’loranes, Christensen-Dalsgaard & Thomp- 

erage seismic parameters (Azz) and icmax- From this sam- 

son 

2005 

Silva Aguirre et al.||2011a 

1 , and are insensitive to 


pie, 35 stars have observations with high enough signal-to- 
noise-ratios (SNRs) for reliable individual frequencies to be 
extracted (see companion paper by Davies et al. 2015, in 
preparation). Typically, more than six consecutive modes 
are detected for each angular degree ^ = 0 — 2 allowing the 
construction of frequency combinations that aim at isolat¬ 
ing the contribution of the inner structure from that of the 
outer layers. In particular, the frequency ratios ro 2 , roi, and 
rio defined as (Roxburgh fc Vorontsov|2003 1 


ro2(n) = 


roi(n) = 


do2{n) 

doi{n) 


Av\{n)' Avo{n+l)' 

where do 2 {n) = Vnfi — yn- 1,2 is the small frequency sepa¬ 
ration and do\{n) and dio(n) are the smooth 5-point small 
frequency separations: 

doi(n) = -I- 6 ^'n ,0 - 4l/n,l + t2n+i,o) (4) 

O 

dio(n) = - 4iz„,o -b - 4i/„+i,o -b I'n-t-qi) .(5) 

O 

The methods considered in this paper to determine 
stellar properties use different sets of asteroseismic observ¬ 
ables in their optimisation procedure (for details see sec¬ 
tions 3.2 and 4.3 below, as well as Appendix [aI) : only the 


dio(r 


( 2 ) 


( 3 ) 


individual oscillation frequencies, only the frequency ratios, 
or a weighted average of both sets. For reference, typical 
uncertainties in the three closest frequencies to iZmax are 


This allows direct comparison of theoretical frequencies with 
the observed quantities without resorting to a surface cor¬ 
rection, such as the power-law term applied to a solar model 


to reproduce the helioseismic observations (e.g., Kjeldsen, 


Bedding & Christensen-Dalsgaard 20081. Moreover, it has 


been shown that use of the frequency ratios provides the 
most precise results on asteroseismic ages as compared to 
htting individual frequencies or using the scaling relations 
( Silva Aguirre et al.|2013 Lebreton &: Goupil||2014 ). 

The ratios and their corresponding covariance matri¬ 
ces have been obtained directly from the MCMC chains 
following the procedure described in Davies et al. (2015, 
in preparation). It is customary to write the ratios roi(n) 
and rio(n) as one unique set of observables, called roio (see, 
Silva Aguirre et al.|2013 |: 


roio = {roi(u), rio(n), roi(n -b 1), rio(n -b 1),...} . 


( 6 ) 


While the frequency ratios provide a reliable proxy of 
the conditions in the deep stellar interior, their formulation 
is based on the assumption that the oscillation frequencies 
behave purely as pressure modes ([Roxburgh fc Vorontsov] 


2000 2003 and references therein). Once a star exhausts 


its central hydrogen supply at the end of the main-sequence 
phase, its core contracts rapidly resulting in coupling be¬ 
tween the propagation cavities of purely acoustic p modes 
and gravity g modes. Consequently, an oscillation mode can 
take mixed character and behave as a g mode in the cen¬ 
tral parts of the star while retaining its p-mode character in 
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Figure 2. Sample distribution in the Tejf v/s Av diagram, de¬ 
picted over stellar evolutionary tracks at solar metallicity for 
masses between 0.8-1.3 Mq. Observational uncertainties in Av 
are smaller than the symbol size. 


the outer layers; this feature is clearly visible in the echelle 
diagram (Grec, Fossat & Pomerantz 19831 as an avoided 
crossing ( Aizenman, Smeyers fc Weigert|1977 1. 

The surface amplitudes of these mixed modes are 
smaller than those of pure p-modes and they are therefore 
harder to detect. Two of the 35 targets show clear evidence 
of modes with mixed character in their i=\ frequencies, for 
which the ratios are not applicable. Their analysis will be 
presented elsewhere using tools especially designed for deal¬ 


ing with these types of oscillation modes (e.g., Deheuvels fc 
Michel|[2011 Benomar et al.|[2012|. A few of the remaining 


33 stars show some hint of mixed modes starting to develop 
or seemingly present in their 1 = 2 ridges. In these cases only 
the quantity roio was matched as described in Section [3.2| 
Figure shows the echelle diagram of two KOIs from our 
sample. The mixed-mode nature in the i = 2 ridge around 
~1070 pHz is very clear in the case of Kepler-36, while the 
ridge is completely smooth (and thus p-mode like) in the 
less evolved star Kepler-409. 

The hnal sample of stars included in this study is shown 
in Fig. together with theoretical stellar evolution tracks. 
Our targets span a wide range of effective temperatures 
(5050 < Teft < 6450 K), metallicities (—0.4 < [Fe/H] < 
-1-0.30 dex), and evolutionary stages that will require the in¬ 
clusion of different input physics to correctly determine their 
intrinsic properties. 


3 STELLAR MODELS AND STATISTICAL 
ANALYSIS 

The following sections describe the physical processes in¬ 
cluded in our stellar models, the Bayesian approach, assump¬ 
tions and observables used to determine the central values 
and statistical uncertainties, and the adopted set of spec¬ 
troscopic constraints of our targets. To test the effects in 
our derived stellar properties of different input physics, as 
well as properly accounting for systematic uncertainties in 


our results, we computed several grids of evolutionary tracks 
changing one physical ingredient at the time. 


3.1 Grids of evolutionary models 


We have computed grids of stellar models using the Garching 


Stellar Evolution Code (GARSTEC, Weiss & Schlattl 20081, 
summarised in Table The basic input physics common to 
all grids includes the 2005 version of the OPAL Equation 
of State (Rogers, Swenson fc Iglesi^llOOO Rogers & Nay- 


fonov 20021, OPAL opacities for high temperatures (Igle- 


sias & Rogers 1996 1 and those of Ferguson et al. (20051 


for low temperatures, the NACRE compilation of nuclear 
reaction rates (Angulo et al.| 1999|) including the updated 


^N(p, 7 ) °0 reaction from 


Formicola et al. 


(20041, and the 


mixing-length theory of convection as described in |Kippen-| 
|hahn, Weigert fc Weiss (20121. 

In our computations the quantity [Fe/H] is defined as 
[Fe/H] — log(2su,./Xsur) log(2sur/Ngur) o, where .^sur and 
Aaur correspond to the surface heavy element and hydro¬ 
gen fractions, respectively. This relation between [Fe/H] and 
ZjX is only valid if the element distribution relative to iron 
is the same in the observed star as in the Sun. A mea¬ 
surement of alpha abundances [a/Fe] is available from high- 
resolution spectroscopy for one of our targets (Kepler-444, 
see |Campante et al.||2015 (. Their obtained average value of 
[a/Fe] = -1-0.26 ± 0.07 has been transformed to ZjX using 
the prescription of |Salaris, Chiefli fc Straniero| ( [T993[ ). 

The initial compositions of our evolutionary models 
were determined using a helium-to-metal enrichment law an¬ 
chored to the Big Bang nucleosynthesis primordial values of 


Zq = 0.0 and Yo = 0.248 (Steigman 20101. A linear ht to 


the initial abundances of a solar calibration results in a re¬ 
lation close to AY IAZ = 1.4, which is the one adopted in 
our grids (see e.g., Pietrinferni et aL]|2004 (. The impact on 
stellar parameters of variations in the ratio between helium 
and heavy elements is explored in Section [4.4[ 

The efficiency of convection was determined from 
a standard solar calibration using microscopic diffusion 
( Thoul, Bahcall fc Loeb||l994p and the [Grevesse fc Sauval| 
(19981 solar mixture, resulting in ajuLT = 1.791. Since 


GARSTEC does not include the effects of radiative levi¬ 
tation or turbulent mixing below the base of the convective 
envelope, we have restricted the calculations of the grid with 


diffusion to masses below 1.2 M© (see Turcotte, Richer fc 
|Michaud|1998p . Our diffusion grid comprises masses between 
0.7 M 0 - 1.2 Mq in steps of 0.01 Mq and initial compositions 
of —0.65 < [Fe/H]< -1-0.50 in steps of 0.05 dex. This mass 
range has been extended in our standard grid (i.e., without 
diffusion) to 0.7 Mq-1.8 M©, also in steps of 0.01 M©. Note 
that we have kept the value of omut constant in these grids 
because solar models calibrated with and without diffusion 
show variations in this parameter of Aomut ~ 0.12 (see, 
e.g., Ciacio, Degl’innocenti fc Ricci||1997 1. Our intention is 
to explore systematic uncertainties by changing one param¬ 
eter at the time, and as part of this effort we explore larger 
changes in aiuLT (see below). 

Furthermore, we constructed a grid including convec¬ 
tive overshooting in the diffusive approach implemented in 
GARSTEC (see Section 3.1.5 in Weiss fc Schlattl|2008 1 , us¬ 
ing an efficiency of ^ = 0.016 that has been calibrated to 


open clusters (e.g.. Magic et al. 2010 Silva Aguirre et al. 
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Table 1. Description of the input physics considered in each of the grids of models. For all grids, mass separation is 0.01 Mq and 
metallicity separation is 0.05 dex. When diffusion is included this is done using the prescription of|Thoul, Bahcall fc Loeb| | |1994[ |. 
Overshooting is incorporated as an exponential decay of the convective velocities (see, [Weiss &: Schlattl||2008| |. Its efficiency has been 
calibrated with open clusters to ^ = 0.016 ijMagic et al.|20Tol|. 


Name 

Mass (Mq) 

[Fe/H] 

Solar Abundances 

Diffusion 

Overshooting 

«MLT 

AY/AZ 

GS98sta 

0.70; 1.80 

+0.50;-0.65 

Grevesse & Sauval 


1998 

1 No 

No 

1.791 

1.4 

GS98ove 

1.00; 1.80 

+0.50;-0.65 

Grevesse & Sauval 


1998 

1 No 

Yes 

1.791 

1.4 

GS98dif 

0.70; 1.20 

+0.50;-0.65 

Grevesse & Sauval 


1998 

1 Yes 

No 

1.791 

1.4 

GS98al+ 

0.70; 1.80 

+0.50;-0.65 

Grevesse & Sauval 


1998 

1 No 

No 

2.091 

1.4 

GS98al- 

0.70; 1.80 

+0.50;-0.65 

Grevesse & Sauval 


1998 

1 No 

No 

1.491 

1.4 


AS09 


0.70 


1.80 +0.50;-0.65 


Asplund et al. (2009 \ 


No 


No 


1.865 


1.4 


2011a I and is roughly equivalent to the aov = 0.2 pressure 


scale heights commonly obtained when instantaneous mix¬ 
ing in the overshoot region is used. Additional grids of mod¬ 
els were constructed changing the value of the mixing-length 
parameter and the reference set of solar abundances (called 
GS98al+, GS98al—, and AS09 in Table[^ see also Section [4^ 
below). The input physics of the grids is detailed in Table 111 
and we will use the naming convention given there for the 
remainder of this paper. 

For hundreds of models along each evolutionary 
track, we calculated theoretical oscillations frequencies us¬ 
ing the Aarhus Adiabatic Oscillation Package (ADIPLS, 


Ghristensen-Dalsgaard 2008aI. Each grid comprises over half 


a million models fully covering the spectroscopic and aster- 
oseismic parameter space of our targets. The final stellar 
properties for our targets have been obtained with either 
the GS98dif or the GS98ove grids, for reasons thoroughly 
discussed in Sections 14.1.11 and 14.1.21 below. 


3.2 Bayesian approach 


The stellar properties of each star are determined using 
an adapted version of the Bayesian approach described in 


Serenelli et al. (20131 that we have named the BAyesian 


STellar Algorithm (BASTA). If v is a set of model stellar 
parameters (e.g. mass, composition, age), and O the ob¬ 
served data, then according to Bayes’ rule the probability 
density function (PDF) of v given O is 


p(v|C>) ocp(v)£(C)|v). (7) 

Here, C (Cl|v) is the likelihood of O given v and p (v) is the 
prior PDF of v representing any prior knowledge on these 
quantities. 

To determine p (v) we considered the following. The 
studied sample is comprised of stars presenting transit-like 
signals in their light curves, as well as stochastically excited 
oscillations. Since the selection function of pulsating exo¬ 
planet candidate host star observed by Kepler is complex, 
we cannot determine a prior probability for stellar proper¬ 
ties based on these characteristics. Therefore, we assume a 
flat prior in [Fe/H] and age including only a strict cut on 
the latter at 15 Gyr (to properly construct the age PDF 
of old stars), and we use a standard Salpeter Initial Mass 
Function. 

We compute the likelihood of the observed values O 
given a set of model parameters v assuming Gaussian dis¬ 


tributed errors. BASTA is flexible in its input, and the ob¬ 
servables included can be the spectroscopic parameters (Teff, 
[Fe/H], log g), the individual oscillation frequencies, the 
global average asteroseismic parameters (Av) and t'max, the 
frequency ratios ro 2 and roio defined in Equations[^and[^ or 
any combination of the aforementioned quantities. The con¬ 
struction of the ratios introduces correlations as a function 
of frequency which are taken into account in the likelihood 
calculation of each model: 


£(0|v) = 


(2^)i/2^ 


exp (-xV2) 


( 8 ) 


^ — (Oobs Omodel) C (Oobs Omodel) , (9) 


where C is the covariance matrix of the observed values 
Ooba, and Oniodei are the same quantities determined from the 
model. For computational efficiency we preselect models for 
determining the likelihood within 500 K and 0.5 dex of the 
spectroscopic Teff and [Fe/H] values, and having an average 
large frequency separation (Av) within 15% of the observed 
one (approximately S-cr from the median uncertainty for the 
majority of the asteroseismic Kepler solar-type sample, see 


Ghaplin et al. (2014l). The value of {Av} in the models is 
obtained using theoretical frequencies of oscillation following 
the prescription of jWhite et al.| | f2011| ) (see section]^ below). 

Our final values for stellar properties of all targets were 
determined with BASTA and GARSTEG grids using Tefr, 
[Fe/H], and the frequency ratios roio and ro 2 as the observ¬ 
ables included in the likelihood calculation (except in a few 
cases where only roio was fitted because £ = 2 mixed modes 
prevented us from using ro 2 , see section . For compari¬ 
son we have also determined all stellar properties using the 
individual oscillation frequencies (see Section |4.3| ) and the 
global average parameters {Aty) and r'max (see Section]^ as 
the asteroseismic data to be fitted. 


Recently, Roxburgh & Vorontsov (20131 proposed that 


fitting the separation ratios roio and ro 2 as a function of 
radial order could lead to biases in the determined stellar 
structure, and recommended comparing observations with 
models by interpolating the theoretical ratios at the ob¬ 
served central frequencies. We have implemented this fit in 
BASTA and found no significant difference with the results 
obtained when comparing at a fixed radial order. Similar 
conclusions were given by [Lebreton &: Goupilj ( |2014| ) (see 
their Appendix A.l). 
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Once the probability density function has been com¬ 
puted using the definitions above, the marginalised posterior 
PDF of any stellar quantity x is obtained as 

p(x\0) — J S{x(y) — x)p{v\0)wv<Av, (10) 

where Wv are the appropriate weights used to account for 
the volume of parameter space occupied by a stellar track 
point characterised by v, i.e. half the distance to its neigh¬ 
bouring points in mass, metallicity, and age. The final values 
reported for each stellar quantity are obtained from the me¬ 
dian of the posterior PDF and the 16 and 84 percentiles. 


3.3 Spectroscopic input parameters 


We initially adopted the spectroscopic values reported in 


Huber et al. (2013b I as input parameters in the Bayesian 


procedure described above. Since this approach allows de¬ 
termination of the posterior PDF from models of any stellar 
quantity, comparison showed that in two stars the best as- 
teroseismic solutions favoured atmospheric parameters not 
compatible with the spectroscopic ones. To explore the rea¬ 
sons for this discrepancy, we have made a new determina¬ 
tion of Teff and [Fe/H] from high-resolution spectra on these 
single targets as available on the website of the Kepler Com¬ 
munity Follow-up Observing Program (CFOFQ. 

The method takes advantage of the accurate values of 
the stellar surface gravity, log g, determined from seismic 
data, which allows a determination of Tefif from the require¬ 
ment that the Fel and Fell lines should give the same Fe 
abundance. Given that [Fe/H] derived from Fel lines in¬ 
creases by ~0.055 dex for solar-type stars when Tefi is in¬ 
creased by 100 K, and that [Fe/H] based on Fell lines de¬ 
creases by ~0.03 dex for the same change in Tefi, the method 
leads to precise values of Teft (±70 K) if the two [Fe/H] val¬ 
ues are determined with errors less than ±0.03 dex. 

This is possible when high-resolution spectra having 
S/N > 40 are analysed differentially with respect to the Sun. 
Based on the solar spectrum we have selected 39 Fel and 14 
Fell unblended lines in the wavelength region 480—650 nm. 
Their equivalent widths (EWs), which range from 30 to 80 
mA in the stellar spectra, were measured by Gaussian fitting 
to the line profiles and analysed with MARCS model atmo¬ 
spheres ( Gustafsson et al.|2008 1 to provide LTE abundances 
as described in Nissen et al. (20141. The Fe abundances from 


Fel lines are subject to small non-LTE effects, which were 
taken into account using the calculations of |Lind, Berge-| 
mann & Asplund (20121. These corrections decrease the de¬ 


rived Teff by 15 to 30 K. 

The abundances are based on medium to strong spec¬ 
tral lines that are sensitive to the microturbulence veloc¬ 
ity Cturb- This parameter is not well constrained from the 
set of lines used, and we have therefore adopted a relation 
^turb (km/s)= 1.0 ±0.00063(reff- 5777) - 0.54(log g-4.44) 
derived for solar-type stars with high S/N HARPS spectra 
( Nissen et al.|20i4 K The derived Tea is not very sensitive to 
^turb, because the two sets of Fel and Fell lines have about 
the same average EW. For an estimated error of 0.2 km/s 
in ^turb, the error in Tea is 25 K, whereas the corresponding 


^ https://cfop.ipac.caltech.edu/home/ 


error in the derived [Fe/H] is 0.05 dex. Furthermore, we note 
that an error of 0.05 dex in log g induces an error of 20 K 
in Teft, and has a negligible effect (<0.01 dex) on [Fe/H]. 

The results obtained are given in Table The listed 
errors correspond to the uncertainties in [Fe/H] estimated 
from the line-to-line scatter of the derived abundances. The 
small errors arising from the uncertainties in log g and ^turb 
were added in quadrature. For the two stars in Table 
Tefi and [Fe/H] deviate significantly from the values deter¬ 
mined in previous studies, i.e. KOI-5 {Tea = 5753 ± 75 K, 
[Fe/H] = 0.05 ±0.1 dex, [Huber et al.|2013b| ), and Kepler-21 


(Tefi = 6131±44 K, [Fe/H] = —0.15±0.06 dex, Howell et al. 


20121. The final stellar properties in this paper were deter¬ 


mined using these updated results. We emphasise that, for 
completeness, we analysed another five stars from thejHuberj 


et al. (2013b I sample and found good agreement: average 


differences and standard deviations for these targets (new 
values— [Huber et al. (2013b l) are ATefi = —13 ± 70 K and 
A[Fe/H] = -0.06 ± 0.06 dex. 


4 RESULTS 


To determine the final stellar properties and uncertainties 
of our sample, we have devised the following procedure. The 
central values and statistical uncertainties were obtained us¬ 
ing the BASTA and a choice of GARSTEG grids including 
different input physics as described in Section [4.1[ These 
are the recommended stellar properties of the sample and 
are listed in Table We then explored the impact of chang¬ 
ing the reference solar abundances and mixing-length pa¬ 
rameter (section [4.2| to estimate a systematic uncertainty 
from the input physics. We then compared our results with 
those obtained from three different asteroseismic pipelines, 
providing a systematic uncertainty arising from the use of 
different evolution codes and fitting algorithms (section [4.3[ ). 
Finally, we estimated a systematic uncertainty arising from 


our choice of AK/AZ = 1.4 (section 4.41, and discuss the 
possibility of breaking the degeneracy between helium and 
mass using Gaia results. 


4.1 Adopted physics for stellar properties and 
statistical uncertainties 

The BASTA described in Section [3.2[ was initially applied 
to the sample using the standard grid (GS98sta, see Ta¬ 
ble]^. Figurej^shows the obtained distribution of masses for 
the targets, encompassing a range between ~0.75—1.6 Mq. 
There are several physical processes that can considerably 
affect the main-sequence evolution of stars, but whose effi¬ 
ciency and range of applicability is not fully constrained by 
theory or independent observations. We explore in the fol¬ 
lowing sections the cases of microscopic diffusion and over¬ 
shooting. 


4-1.1 Microscopic diffusion 

In the low-mass regime, neglecting the effects of micro¬ 
scopic diffusion can have significant impact on the deter¬ 
mined properties. If stars are indeed subject to element set¬ 
tling throughout their main-sequence lifetimes, the derived 
ages can be particularly biased (e.g., Valle et al. 2014b). 
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Table 2. Spectroscopic determinations of selected KOIs. The log g values were adopted from asteroseismic determinations of |Huber| 
|et al.| |2013b[|. Uncertainties quoted in and [Fe/H] correspond to the internal precision of the method. 


Name 

Spectrum 

S/N 

Teff (K) 

log g (dex) 

[Fe/H] (dex) 

^turb (km/s) 

KOI-5 

McD2.7m 

70 

5945±60 

4.05 

+0.17±0.05 

1.3 

Kepler-21 

TRES 

150 

6305±50 

4.02 

-0.03±0.05 

1.6 



Mass Distribution (M^) 


Figure 3. Mass distribution of the sample obtained using the 
standard grid (GS98sta). 


Evidence for the occurrence of microscopic diffusion is clear 
for the Sun (e.g., Guzik fc Cox|1993 Christensen-Dalsgaard, 


Profhtt fc Thompson|jl9^ l and has been growing in recent 

years for star clusters (see e.g., jKorn__e^_ay [2007| |Gruyteri| 


et al.|2013 Onehag, Gustafsson fc Korn|2014 1 revealing the 

importance of this process in low-mass stars. 

Using the grid of models that includes the effects of dif¬ 
fusion (GS98dif) we have redetermined the stellar properties 
of all stars with masses predicted from the standard grid to 
be below ~1.2 M©. The top panel of Fig.|^shows the result¬ 
ing fractional mass differences, with the two determinations 
typically agreeing within ~2.5%. There is no visible trend 
as a function of metallicity although the largest differences 
are found for the most metal-poor stars, as expected (see 
e.g., Valle et al.|2014b l. However, the ages predicted by the 
standard grid are systematically older than those based on 
diffusion (bottom panel in Fig. |^. Since the masses deter¬ 
mined from both grids agree on average, this behaviour is 
not what would be expected from a pure mass effect and 
is a consequence of several physical effects combined that 
balance each other. 

Stellar models that include microscopic diffusion sink 
helium and heavy elements while diffusing hydrogen towards 
the surface during their main-sequence lifetime. This process 
effectively reduces the amount of fuel available in the core 
producing a faster main-sequence evolution and therefore 
decreasing the age a model has at a given mean density. 
Interestingly the effects of diffusion on age counterbalance 
those expected from a pure mass effect, since even if the mass 
predicted by the non-diffusive model is higher than the dif¬ 
fusive one its age remains older. Thus, despite not producing 
biases in the masses determined, the use of diffusion does af- 



[Fe/H] 


Figure 4. Fractional differences in mass and age obtained from 
the standard and diffusion grids of models as a function of spec¬ 
troscopic metallicity, in the sense (GS98sta-GS98dif) Top: mass 
comparison. Bottom: age comparison. See text for details. 


feet the ages obtained from asteroseismic fitting and should 
be taken into account in the relevant regime. Similar results 


were found by e.g., Miglio & Montalban (20051; Valle et al. 
| |2014a| ), where models with diffusion predict younger ages 
than those not including this effect. 

Gompared to canonical models, the transport of chem¬ 
ical species by microscopic diffusion produces different 
composition profiles and therefore different opacity values 
through the stellar interior. As a consequence, evolution¬ 
ary tracks computed with the same initial chemical mixture 
including or not these effects will also show a different effec¬ 
tive temperature evolution during the main-sequence phase. 
An example of the effects played by microscopic diffusion 
is shown in Fig. These results are obtained for Kepler-93 
using the frequency ratios roio, ro 2 , and spectroscopic pa- 
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Figure 5. Values of determined using the frequency ratios and spectroscopic constraints for Kepler-93. Solid lines depict the 1 a 
uncertainties in the spectroscopic atmospheric parameters. Left standard grid. Right, diffusion grid. 


rameters as the observational data to be reproduced. The 
plots depict the surface from each grid of models in 
the Tefi versus [Fe/H] plane as well as the observed spec¬ 
troscopic properties and uncertainties. Due to the number 
of observables and size of the fractional uncertainties, the 
fit is dominated by the asteroseismic data and the grid ne¬ 
glecting diffusion favours a parameter space that is not fully 
compatible with the observed Tefr and [Fe/H] of the target. 
In particular, GS98sta results predict Teff values too high 
compared with the spectroscopic ones. For this reason, the 
grid with diffusion has been used for determining the final 
central values and statistical uncertainties of targets with 
masses below ~1.2 Mq. 

One must be cautious when determining uncertainties 
using the diffusion grid because it suffers from an unavoid¬ 
able edge effect for masses above ~1.15 Mq, where the pos¬ 
terior PDF are artificially cut due to the chosen upper mass 
limit when constructing the grid (1.20 Mq). To ensure a 
smooth transition in the obtained stellar properties when 
switching from GS98dif to another grid, we compared the 
results predicted with GS98sta and GS98ove to those of 
GS98dif in six targets at the edge of this mass range (KOI-5, 
Kepler-65, Kepler-25, Kepler-126, KOT268, and Kepler-129, 
see Table 1^. All stars where the GS98dif grid returned cen¬ 
tral mass values above 1.15 Mq clearly showed the edge 
effect in the uncertainties, but the central values in the stel¬ 
lar properties are fully compatible with those obtained from 
the other two grids (note that the effects of overshooting are 
not expected to be significant in the 1.15-1.2 Mq range, see 
Section 4.1.2 below). The reason is that for these masses, 
the evolutionary timescale is short enough for the effects of 
diffusion to be minor. Thus, in our recommended values pre¬ 
sented in Table we restrict ourselves to the diffusion grid 
in all stars with central mass values below 1.15 M©. 


4-1.2 Convective overshooting 

In stars with masses higher than ~1.1 Mq the size of the 
convective core severely affects the age at which central hy¬ 
drogen exhaustion (the terminal age main-sequence, TAMS) 
is reached. The extent of the central mixed region is con¬ 


trolled by the convective overshoot efficiency, a free param¬ 
eter in stellar evolution calculations. There is a wealth of 
evidence in the literature that convective core sizes are prob¬ 
ably larger than those predicted by models applying the pure 
Schwarzschild criterion for convective boundary determina- 


tion (e.g., Maeder|1974b|a 

Maeder & Meynet|1991 Ghiosi, 

IBertelli & Bressan||1992t 

Vandenberg et al. 1120071 .[ust to 


name a few), and recent analyses of AepZer stars have further 
supported this scenario (Silva Aguirre et al. [20131. Thus, we 


have applied our Bayesian method using the overshooting 
grid (GS98ove) for all targets and compared the results to 
those from the standard set of models. 


Figure shows the fractional difference in mass and 
age obtained from both grids for stars with masses above 
1.1 Mq. In three targets (Kepler-145, KOI-974 and Kepler- 
21 ) the mass and age deviations can be attributed to the 
quality of the fits: the PDFs of these particular stars show a 
clearly bimodal distribution (see below). A similar situation 
is observed in the target at ~1.5 Mq, where the large mass 
and age uncertainty in the GS98sta results is a consequence 
of its Bayesian results encompassing the secondary peak in 
the distribution. 


Besides the obvious outliers (red squares in Fig. §, dif¬ 
ferences in the results are noticeable in five targets (HAT- 
P7, KOI-5, Kepler-50, Kepler-129, and KOI-288). Figure 
shows the fractional mass difference as a function of the 
TAMS ratio, defined as the ratio between the age of each 
target’s best fit model and the age that model has when 
exhausting its central hydrogen content. The differences in 
mass appear at TAMS ratio values of ~0.9 when convec¬ 
tive cores are well developed and strongly influence the fre¬ 
quency ratios (Silva Aguirre et al. 2011al. Moreover, once 
the stars reach the turn-off at TAMS ratio of ~1.0 some 
results are fully compatible while others still show small dif¬ 
ferences. This might be an indication of asteroseismic data 
keeping record of past mixing processes in stars once these 
have evolved slightly beyond the main-sequence phase (see 
Deheuvels fc MichelpOlf where similar results were found). 


We note that for our targets with TAMS ratio larger than 
1 , their best-fit models show evidence of mixed modes with 
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1.1 1.2 1.3 1.4 1.5 1.6 

Mass {Mq) 


Figure 6. Fractional differences in mass and age obtained with 
the standard and overshooting grids as a function of stellar mass, 
in the sense (GS98sta-GS98ove) Top: mass comparison. Bot¬ 
tom: age comparison. Red squares shows the positions of outliers 
Kepler-145, KOI-974, and Kepler-21. See text for details. 


no obvious counterparts in the observations, probably due 
to their low amplitudes. 

The three outliers present in Figs. and (Kepler-145, 
KOI-974, and Kepler-21) stem from bimodal probability dis¬ 
tribution functions obtained from the GS98sta grid, which 
are not observed in the results when the GS98ove models 
are used instead. Figure shows the example of KOI-974, 
where the GS98sta grid favours a star with mass ~1.27 Mq 
that has evolved beyond the main sequence (TAMS~1.1, 
see Fig. [^. However, there is evidence of a secondary solu¬ 
tion at higher masses, close to ~1.43 Mq. When using the 
overshooting grid instead the PDF shows a unique peak at 
~ 1.39 Mq compatible with the high-mass solution seen in 
the GS98sta grid. A similar behaviour is observed in the 
PDFs for Kepler-145 and Kepler-21, where the GS98ove 
results favour masses consistent with the higher-mass sec¬ 
ondary peak in the GS98sta probability distribution func¬ 
tion. We note that the bimodal nature of the solution for 
Kepler-21 had already been pointed out by |Howell et al.| 
1 2 ^ . 

For the three outliers the inclusion of overshooting pro¬ 
duces a better ht to the frequency combinations and argues 
against the lower-mass peak in the PDF. Since the GS98sta 



Figure 7. Fractional mass difference between the standard and 
overshooting grids as a function of the Terminal Age Main Se¬ 
quence (TAMS) ratio for the GS98sta grid. Red squares shows 
the positions of outliers Kepler-145, KOI-974, and Kepler-21. The 
TAMS ratio is defined as the age of the best fit model in the 
GS98sta grid divided by the age at the end of the core hydrogen¬ 
burning phase. See text for details. 



o.oi—- 

1.15 1.20 1.25 1.30 1.35 1.40 1.45 1.50 1.55 

Mass (Mg) 


Figure 8. Mass distribution obtained with BASTA for KOI-974 
using the standard and overshooting grids. 


and GS98ove grids cover the same parameter space in terms 
of initial masses and compositions, we also compared the 
average likelihood of the grids for these targets and always 
found better agreement with the overshooting models. In 
fact, the GS98ove results in all three cases favours stars 
close to the end of the main-sequence phase (TAMS ~0.92), 
where convective cores are well developed and masses de¬ 
termined with the standard grid are slightly higher than 
with overshooting (see Fig. [^. Thus, the final parameters 
of these three KOIs, given in Table are those determined 
from the GS98ove. This solution for Kepler-21 also results 
in a luminosity and distance compatible with that deter¬ 
mined from parallaxes, as described in section |4.4| below. 
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Fractional Uncertainty (%) 


Figure 9. Fractional uncertainties added in quadrature in the 
final stellar parameters determined with our BASTA. See text 
for details. 


The actual mass value at which extending the convective 
core size by means of overshooting has an impact in the ob¬ 
tained properties is expected to be between 1.2 — l.SM©, 
but the quality of the data and mass coverage of the current 
sample are not adequate to thoroughly test this transition. 
Figures and show that this reasonable expectation is 
consistent with the data. 

Based of the analysis and comparisons on the input 
physics discussed in Sections |4.1.1| and |4.1.2[ we choose the 
GS98dif grid to determine the final stellar properties and 
statistical uncertainties for stars of masses below ~1.15 M©, 
and the GS98ove grid for the remainder of the sample. These 
values are given in Table In terms of internal precision, 
the median fractional uncertainties returned by the BASTA 
are 1.2% (radius), 1.7% (density), 3.3% (mass), and 14% 
(age). We note that the seemingly large fractional uncer¬ 
tainty on the density compared with asteroseismic determi¬ 
nations made by other methods arises from the use of fre¬ 
quency ratios instead of individual oscillation frequencies. 
The distributions for our sample shown in Fig. [^correspond 
to the statistical uncertainties only. To check how robust this 
error budget is we explore in the next sections the system¬ 
atic contribution to the uncertainties from the input physics 
and fitting algorithms. 


realistic than the simple ID atmospheres used by |Grevesse| 


& Sauval (19981, their results are not in agreement with 


the helioseismic determinations of the solar internal struc¬ 
ture (e.g., Basu fc Antia|2(i08 Christensen-Dalsgaard et al. 


2009 Serenelli et al.|2009 1. Similarly, these simulations pre¬ 

dict variations in the efficiency of convection as a function of 
Tefi and [Fe/H] (e.g., Trampedach et al.|2014 Magic, Weiss 


fc Asplund||2014|). In hydrostatic evolutionary models this 


process is mimicked by the mixing-length theory and its ef¬ 
ficiency is represented by one parameter (omlt) calibrated 
to reproduce the properties of the Sun. 

A thorough analysis of the virtues of 3D model pre¬ 
dictions compared to those of ID models goes beyond the 
scope of this paper, but we are interested in knowing the 
impact of these results on our fitted stellar properties. For 
the case of the solar abundances we have chosen to use the 


compilation of Asplund et al. (20091 since its determinations 


represent an extreme case with respect to those of|Grevesse| 
|fc Sauval ( |1998[ ), while the results of |CafFau et ah (20111 
lie approximately in between these two determinations. In 
terms of convective efhciency our standard grid uses a value 
calibrated to the solar properties; to test the impact of only 
changing omlt we have computed grids with the|Grevesse fc| 


Sauval 119981 abundances varying this parameter by ±0.3, 


in agreement with the changes predicted by 3D simulations 


(see e.g., Trampedach et al. 20141. The resulting sets of mod¬ 
els (AS09, GS98al±, GS98al—) are summarised in Table[^ 
The stellar parameters obtained with BASTA and these 
three grids deviate systematically from the GS98sta re¬ 
sults. When comparing to the values obtained with the |As-| 
[plund et al.| ( |2009| solar abundances, the median and stan¬ 
dard deviation of the fractional differences is —0.1 ± 0.3% 
(density), ±0.3 ± 0.3% (radius), ±0.9 ± 0.6% (mass), and 
±4.0 ± 3.3% (age). An increase in the mixing-length pa¬ 
rameter on the other hand produces changes at the level of 
±0.9 ±0.7% (density), -0.9 ± 0.6% (radius), -2.4 ± 2.2% 
(mass), and—0.2 ± 9.0% (age), and of opposite sign if the 
Q!mlt value is decreased. Given that these differences are 
always smaller than the statistical errors quoted in the pre¬ 
vious section, we have not added them to the error budget 
given in Table [^ We consider the dispersions as the major 
contribution to the uncertainties, and add them in quadra¬ 
ture to quote a reference systematic effect arising from the 
difference in the input physics of 0.8% (density), 0.7% (ra¬ 
dius), 2.3% (mass), and 9.6%(age). 


4.2 Systematic uncertainties from input physics 


Our central recommended values and statistical uncertain¬ 
ties derived in the previous sections were based on fixed sets 
of input physics chosen according to our current best knowl¬ 
edge of stellar evolution (i.e., including microscopic diffusion 
and core overshooting in the relevant mass regimes). In this 
section we explore the impact on our results from changing 
two quantities whose values are not yet fully constrained by 
observations or theory, namely the solar abundances and the 
mixing-length parameter. 

The photospheric composition of the Sun has been re¬ 
vised in recent years inspired by the predictions of 3D model 


atmospheres (e.g., Asplund, Grevesse & Sauval 2005 As- 
[plund et al.||200^ [Gaffau et al.||2011| ). Although these hy 
drodynamical models of the solar atmosphere are much more 


4.3 Systematic uncertainties from fitting methods 

As described in Sections |3.2| and |4.1[ our recommended stel¬ 
lar properties are derived fitting the frequency ratios roio 
and ro 2 and the spectroscopic parameters to two chosen 
grids of models depending on the mass range. In many aster¬ 
oseismic analyses it is customary to fit the individual modes 
frequencies (after applying a suitable surface correction) in¬ 
stead of frequency combinations. To isolate the impact of 
the choice of fitted seismic observables, we take advantage 
of the flexibility of BASTA in its input and compute stellar 
properties fitting the spectroscopic Teg and [Fe/H] and the 


^ Giving the scatter from changing the mixing-length half the 
weight as its contribution is symmetric 
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individual oscillation frequencies (corrected with the Kj eld- 
sen, Bedding & Christensen-Dalsgaard (|2008|) prescription). 


The median fractional differences in the results obtained 
with BASTA+GS98sta using frequencies and ratios is below 
0.5% in density and radius, 1% in mass and 7% in age. These 
values should be considered as a minimum set of systematic 
uncertainties based only on the choice of observables. 

Applying the BASTA to one grid of GARSTEC mod¬ 
els corresponds to a single determination of stellar proper¬ 
ties at a Hxed set of input physics, evolutionary code, pul¬ 
sation code, and minimisation technique. Throughout the 
years several algorithms have been developed to reproduce 
asteroseismic data making use of a variety of stellar evolu¬ 
tion and oscillation codes, considering different assumptions 
on the physics included and the treatment of the data to 
be fitted. We aim now at determining the systematic un¬ 
certainty arising from the use of different fitting pipelines, 
including the numerics of the evolution and pulsation codes 
as well as the minimisation routine and chosen data sets. 

Bearing this in mind, we have obtained results for our 
sample of KOIs using three different methods: the ASTEC 
Fitting (ASTFIT) method, the YALE Monte-Carlo Method 
(YMCM), and the Asteroseismic Modeling Portal (AMP). 
To isolate the effects of using different codes and assump¬ 
tions on data treatment each algorithm has been applied 
using input physics as similar as possible to the standard 
GARSTEC grid (GS98sta), i.e., not including the effects of 
microscopic diffusion nor overshooting. Thus, all compar¬ 
isons in this Section are made between the results of these 
htting algorithms and the BASTA-|-GS98sta combination. 
A description of the methods can be found in Appendix [X} 
but we repeat here the main differences in each of the im¬ 
plementations with respect to the BASTA-|-GS98sta. 

The ASTFIT method uses the ASTEC evolutionary 
code (Christensen-Dalsgaard 2008bI and was applied with 
the same physics as in the GS98sta grid (including a solar- 
calibrated mixing-length and /SYjAZ = 1.4). Its main dif¬ 
ference with the BASTA are the asteroseismic variables to 
be fitted: instead of using frequency ratios, this algorithm 
matches the individual frequencies mimicking the surface 
effects via a scaled function determined from the frequency 
corrections found in the solar case (see Appendix |A1| for 
details). 

The YMCM differs from the BASTA-|-GS98sta combi¬ 
nation in its evolution code (YREC, [Demarque et al!]|2008 l 
and its pulsation code (hrst used in [Antia fc Basu 19941. 
Also, the assumptions on the initial composition are differ¬ 
ent: no helium to metal enrichment ratio is applied and the 
initial helium abundance is a free parameter in the min¬ 
imisation process. The stellar properties are obtained from 
a weighted average of the fit to the individual frequencies 
(with a solar-function correction for surface effects, see Ap¬ 
pendix [A^, the ratios, and the spectroscopic constraints. 
YMCM takes into account the correlations introduced by 
the construction of the frequency ratios (see equations 
and[^. 

In the case of the AMP the evolution and pulsation 
code are the same as in ASTFIT ( [Metcalfe, Creevey &l\ 
Christensen-Dalsgaard|2009 1. However its standard configu¬ 


ration adjusts both the mixing-length parameter and initial 
abundances when searching for the best ht model. The latter 
is found using a weighted average between the fit to the in¬ 


dividual frequencies (corrected with the Kjeldsen, Bedding 
fc Christensen-Dalsgaard] ([2008 1 formulation), the frequency 


ratios roio and ro 2 and the spectroscopic constraints (see Ap¬ 
pendix A3l. Correlations are not taken into account in the 
frequency combinations and these are assumed to be inde¬ 
pendent. 

Figure shows the fractional differences obtained in 
mean density and radius from the different methods. It is 
worth noticing that the YMCM could not provide a solu¬ 
tion in four cases where the best fit was found for helium 
abundances below the standard Big Bang nucleosynthesis 
(SBBN) value of Yb = 0.248. There is overall good agree¬ 
ment in the results, with the AMP values showing a larger 
scatter than the other two methods. The closest matches are 
obtained between the BASTA and ASTFIT; not surprisingly 
since these are computed with the exact same input physics. 
Median differences and standard deviations in densities and 
radii are below the ~1% level, while the comparison with the 
YMCM and AMP leads to median differences of the same 
order but a larger scatter (~2%). This is a remarkable level 
of agreement considering the variety of evolutionary and pul¬ 
sation codes used as well as the treatment of asteroseismic 
data. 

There is evidence of a metallicity dependence in the den¬ 
sity and radius differences between the BASTA and YMCM 
results, while the comparison with the AMP does not sug¬ 
gest any correlation with the spectroscopic parameters. Sim¬ 
ilar features can be seen in Fig. HD where the fractional 
differences in mass and age are shown. Once again there is 
excellent agreement between the BASTA and ASTFIT, with 
median and standard deviations in mass and age of 0.8±2% 
and 1.2 ± 9% respectively and within the uncertainties in 
practically all cases. The previous trends with spectroscopic 
parameters remain with respect to the YMCM and AMP, 
namely a metallicity dependence in the first case and seem¬ 
ingly random scatter in the second one. This leads to mass 
estimates differing by up to ~15% and, via the mass-age 
correlation, to corresponding age differences of more than 
40%. 


While the origin of these discrepancies will be explored 
in the following section, for the sake of the current analysis 
we use the scatter in the comparison between the BASTA 
and the ASTFIT results as a measurement of the sys¬ 
tematic uncertainty arising purely from evolutionary codes 
and htting methods, which conservatively accounts for 1% 
in density and radius, 2% in mass, and 9% in age. This 
level of agreement is comparable to the one obtained with 
BASTA-|-GS98sta when changing the observables from ra¬ 
tios to individual frequencies, showing that remarkably pre¬ 
cise determinations of stellar properties are possible from 
asteroseismology despite the use of different evolutionary 
codes as long as the considered input physics is kept as sim¬ 
ilar as possible. 


4.4 Systematic uncertainties from the initial 
helium abundance 

To understand now the origin of the discrepancies between 
YMCM, AMP, and the BASTA-|-GS98sta combination high¬ 
lighted in the previous section, we naturally focus on the 
parameters that have been treated differently in the algo¬ 
rithms. The bottom panels in Fig. show the fractional 
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Figure 10. Fractional differences in stellar density and radius for the three fitting methods compared to the BASTA results as a function 
of spectroscopic metallicity, given as (BASTA-other). 


mass difference as a function of the difference in the ini¬ 
tial helium abundance. In both cases a correlation between 
AM/M and AYini can be seen, with higher masses predicted 
by the YMCM and AMP for lower values of the initial he¬ 
lium fraction. A similar trend can be seen when compar¬ 
ing the luminosities (middle panels in Fig. 121, which turns 
out to be a consequence of the different weights each fit¬ 
ting method gives to the spectroscopic constraints (see Ap¬ 
pendix]^ and section 3.2 for details). 


Since there is good agreement in the median density 
and radii determined by all pipelines (better than 2% with 
a larger scatter for AMP), differences in luminosity arise 
from differences in effective temperature of each model. The 
top panels of Fig. show the match to the spectroscopic 
Tefi value predicted by each algorithm, revealing a very close 
agreement between these codes and the spectra (mean differ¬ 
ence of 20 K and 35 K for AMP and YMCM, respectively). 
ASTFIT and the BASTA on the other hand predict Teft val¬ 
ues systematically higher than the spectroscopic ones at a 
level of ~45 K and 70 K, respectively. The reason for this 
difference is that the AMP and YMCM have the flexibility 
to fit Teff at a fixed radius by changing the luminosity via 
its strong dependence on the molecular weight and thus the 


helium abundance. Our final results from the BASTA given 
in Table 1^ do not suffer from this systematic offset in Teff as 
they include microscopic diffusion which effectively compen¬ 
sates for this effect as shown in Fig. (see also discussion in 
section 4.11. 


The initial helium abundance in low-and intermediate- 
mass stars is a poorly constrained quantity as helium lines 
are not detectable in their spectra. Usually these abundances 
are computed assuming a linear enrichment law AY/AZ an¬ 
chored to primordial values (Yo, Zq) determined from SBBN 
or a solar calibration. The actual slope of the relation has 
been determined from different sources such as low-main- 
sequence stars in the solar neighbourhood ( [Jimenez ] [2003l 
Casagrande et aT)]|2007 l, galactic and extragalactic H-II re¬ 
gions ( Balser|2006 1 , and stellar evolution theory for a given 
initial mass function (e.g., jChiappini, Renda fc Matteuccij 
20021. There is no consensus yet on the relation being linear 


across all metallicity values or the slope changing for dif¬ 
ferent types of stars and metallicity ranges. However, there 
is broad agreement on the expected range for the relation 
being between 1 ^ AY /AZ ^ 3 (except for some very pecu¬ 
liar objects such as uj Centauri, see |Norris|2004| ). Depending 
on the chosen solar composition, determination of the ini- 
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Figure 11. Fractional differences in stellar mass and age for the three fitting methods compared to the BASTA results as a function of 
spectroscopic metallicity, given as (BASTA-other). 


tial helium abundance of the Sun agrees with a slope of 
1.7 ^ AYjAZ ^ 2.2 (see Serenelli fc Basu||2010 l. 

Figure [T3| compares the predictions of helium to metal 
enrichment laws of different slopes with the initial abun¬ 
dances determined by the YMCM and AMP. We remind 
the reader that the BASTA and ASTFIT results are an¬ 
chored to AY/AZ = 1.4 and therefore follow a straight line 
if plotted in the figure. As expected, a large majority of the 
YMCM results lie below this line and in many cases very 
close to the primordial SBBN value which was used as a 
hard edge in the modelling procedure. The AMP results on 
the other hand only force a hard cut in Yni = 0.22 and show 
a much larger scatter as well as five targets with initial he¬ 
lium abundances below the SBBN line. This occurs for stars 
both at low metallicities and at and above the solar value 
{Zini ~ 0.018). The AMP also predicts high initial helium 
for stars of metallicity lower than the Sun, which occurs only 
for one case in the YMCM results. 

Only two of our targets have available Hipparcos par¬ 


allaxes from Van Leeuwen (20071: Kepler-21 (tt = 8.86 ± 
0.58 mas) and KOI-3158 (tt = 28.03 ± 0.82 mas). Using the 
asteroseismic log g value and spectroscopic metallicity, we 


applied the InfraRed Flux Method (see section 4.51 to de¬ 


termine the bolometric flux and the luminosity of the stars. 
The resulting values of L/Lq = 5.54 ±0.78 (Kepler-21) and 
L/Lq = 0.38 ± 0.03 (KOI-3158) are compatible with those 
determined with the BASTA (see Table[^. ASTFIT predicts 
the correct luminosity for KOT3128 {L/Lq = 0.35 ± 0.02) 
while its determination for Kepler-21 results in a slightly 
lower value than observed {L/Lq — 4.59 ± 0.11). The rea¬ 
son for this discrepancy is that ASTFIT favours the low- 
mass solution in the bimodal mass distribution observed for 
Kepler-21 (see section 4.1 1 ; the final results of the BASTA 
are based on the overshooting grid that favours a higher- 
mass value. For comparison the low-mass solution obtained 
with the BASTA±GS98sta predicts a luminosity of L/Lq = 


4.80 


+0.12 


The results from AMP return similar luminosities as 
ASTFIT: L/Lq = 0.36 ± 0.02 (KOI-3158) and L/Lq = 
4.37 ± 0.19 (Kepler-21). However, in both cases the Fjni 
values from AMP are slightly below the SBBN one: the 
luminosity of KOT3158 is still well reproduced because 
at such low metallicities the initial helium abundance at 
AY/AZ = 1.4 is not very different from the primordial one 
(Yini ~ 0.258). The resulting mass of Kepler-21 is higher 
than that from ASTFIT but the luminosity is lower as a 
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Figure 12. Comparison of stellar properties as a function of the difference in the initial helium abundance, given as (BASTA-code). 
Top panels: difference in effective temperature between the spectroscopically measured one and that extracted from the fitting methods. 
Middle panels: fractional differences in the luminosities predicted by the codes. Bottom panels: Fractional mass differences between codes. 
See text for details. 


consequence of the difference in initial helium abundance at 
compositions closer to solar (Yini ~ 0.270). 


The correlation between mass and helium abundance 
severely affects asteroseismically derived stellar properties 
(see section 4.3 in Lebreton & Goupil 20141, and has re¬ 
currently appeared in studies fitting individual frequencies 
(e.g., [Mathur et al.||201^ [Metcalfe et al.|2014p . Assessing if 
the sub-SBBN results have astrophysical meaning goes be¬ 
yond the scope of this paper, but the dependence on helium 
abundance of asteroseismic solutions needs to be further in¬ 
vestigated and thoroughly understood. In order to estimate 
an uncertainty arising from the changes in initial helium 
content, we compare the results of BASTA, YMCM, and 
AMP for the subset of targets falling within the expected 
range of 1 ^ AYjAZ ^ 3. There are five and nine tar¬ 
gets fulfilling this criterion in the YMCM and AMP results, 
respectively (see Fig. |13[ ). The average standard deviations 
around the median of the results are 1.6% (radius), 1.7% 
(density), 3.6% (mass), and 16.8% (age). These uncertain¬ 
ties are comparable to the statistical ones obtained with 


BASTA-l-GARSTEC grids, as described at the end of sec¬ 
tion |4T] 

In an initial attempt to solving the issue of asteroseismic 
analysis favouring sub-SBBN helium models, we consider 
the correlation of properties with initial helium abundances 
shown in Fig. |12| Our results for the two targets where 
parallax-derived luminosities are available suggest a better 
agreement with grids computed assuming AYjAZ = 1.4, 
although the uncertainties are still rather large and AMP 
results are only slightly more than la away in Kepler-21. 
We expect much tighter observational constraints on the lu¬ 
minosity from the Gaia mission that will essentially trans¬ 
late into a constraint on Yni, breaking the degeneracy and 
subjecting the AYjAZ relation to an additional test from 
asteroseismology. 


4.5 Distances 


We have determined distances to our targets combining the 
results of the InfraRed Flux Method ||Casagrande et al.| 


20101 with asteroseismically determined stellar parameters 
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Figure 13. Initial metallicity versus initial helium abundance for the results determined by two different codes. Horizontal dashed 
lines shows the primordial helium abundance predicted by SBBN, while solid lines depict different galactic enrichment laws anchored to 
^0 = 0.248 and Zq = 0. Solar symbol depicts the initial solar abundance from [Serenelli Sz Basu| | |2010[ |. Left: YMCM. Right: AMP. See 
text for details. 


( Silva Aguirre et al.|20lfb |. This technique has been shown 
to provide accurate ffipparcos-quality distances for individ¬ 
ual targets | |Silva Aguirre et al.|2012| ) and has been recently 
applied to the first Galactic Archaeology study of the Kepler 
held (Casagrande et ^|2014l. 


In its current implementation the method relies on 
multi-band griz photometry from the original Kepler In¬ 
put Catalogue (KIC, Brown et al. pMT| , JHKs magnitudes 
from 2MASS, and the spectroscopically determined metal¬ 
licity to recover the bolometric hux as a function of surface 
gravity. The precise log g values obtained from asteroseis- 
mology are included in an iterative process until convergence 
in reddening and effective temperature is achieved (see sec¬ 
tion 3 in Silva Aguirre et al.||2012 for details). 


Unfortunately only 22 of our 33 targets have reliable 
KIC photometry taken in all griz hlters, and we have initially 
determined distances for this subsample using the combina¬ 
tion of the IRFM and asteroseismology. With the aim of pro¬ 
viding a complete and homogeneous set of stellar distances 
for the full cohort, we have used the J — Ks angular diameter 


calibration of Casagrande, Portinari & Flynn (20061 to ob¬ 


tain distances to all targets, and compared them to the full 
IRFM results. For the 22 targets, the median and standard 
deviation of the fractional distance difference between both 
determinations is 0.3 ± 5.1%, showing that the calibration 
works well for the selected sample. The distances given in 
Table|^for all targets have been determined from the colour- 
angular diameter calibration, and have a fractional median 
uncertainty of 4.4%. 

The two stars in our sample with measured Hippareos 
parallaxes allow for an independent check on our distances. 
The resulting values of d = 112.9 ± 7.4 pc for Kepler-21 and 
d = 35.7 ± 1.04 pc for KOI-3158 are compatible within their 
uncertainties to those determined by our method (114 pc 
and 33 pc respectively, see Table [^. 


5 COMPARISON WITH ASTEROSEISMIC 
SCALING RELATIONS 

The analysis presented in the previous sections has been 
possible for a subsample of almost half the currently stud¬ 
ied KOIs where asteroseismic data are available. In the rest, 
SNRs are too low or the frequency resolution not sufficiently 
high to robustly extract individual oscillation frequencies 


(see e.g., Chaplin et al. 20141. Under these circumstances, 
asteroseismic determination of stellar parameters is carried 
out based on the properties of the readily extractable aver¬ 
age global parameters (Ai/) and i^max. The former quantity 
is a proxy of the mean stellar density p ( Ulrich|1986 1 , while 
the latter has been shown to scale with the surface gravity 


and effective temperature (Brown et al. 

1991 Kjeldsen & 

Bedding|1995 Bedding & Kjeldsen|2003 

I: 


( (A^) V 


^max 

i^max,© 


M 


P0 ’ 
R© 


TeS 

Feff,© 


- 1/2 


( 11 ) 


( 12 ) 


Here, p© and Teu,© correspond to the solar properties. 
Consequently, (A;/) and i^max are the basis of the asteroseis¬ 
mic scaling relations 


M 


((Au) ( n, 

i^max,© ^ / \^(Ai/0)y yTeff,0 


R ^ f i2^ax \ / {AU) \ Y Teff ^ 
R 

© y^max,© J v(Ai^©)y yreff,0_ 


(13) 


(14) 


where (Ai/©) and I'niax,© are the solar values as determined 
by the same method used to analyse the data. 

When a determination of Tgu is available, equa¬ 
tions and provide masses and radii independent of 


any stellar models (e.g., Stello et al. 2008 Silva Aguirre 
|et al.|[2011b[ ). Coupled to evolutionary tracks the astero¬ 
seismic scaling relations can also provide an age estimation 
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via the so-called grid-based method, that matches the spec¬ 
troscopic constraints and measured (Av) and r^max to the 


same values predicted from stellar models (e.g., Basu, Chap- 


lin fc Elsworth||2010[ |Gai et al.]|2011[ |Silva Aguirre et al. 


2014 


Lundkvist, Kjeldsen & Silva Aguirre 20141. These 


techniques have been extensively used in characterisation 
of the CoRoT and Kepler samples (e.g., Mosser et al.poio 


Huber et al.|2014l, studies of open clusters (e.g., Basu et al. 


2011 Stello et al. 20111, galactic archaeology (e.g., Miglio 


et al.|2013b Casagrande et al.|2014 1, and will form the basis 


of asteroseismic analysis for the vast majority of targets ob¬ 


served by the repurposed Kepler mission (K2, Howell et al. 

20141 ). 

has direct theoretical 


11 


The {Av} relation in equation 
foundation while the r'max scaling in equation|12|is mostly an 


empirical relation (although see Belkacem et al.||2011|. For 


these reasons the asteroseismic scaling relations need to be 
thoroughly tested and calibrated, and have been the subject 
of scrutiny in the past few years. Studies have shown that 
radii of main-sequence stars predicted from asteroseismology 
agree with determinations from parallaxes and interferom- 


etry (e.g.. 

North et al. 

2007 

Bedding et al.[[2010 

Miglio 

20111 

Silva Aguirre et al.||2012| |Huber et al.||2012 

White 

et al. 

2013[), while asteroseismic masses still remain to be 


confirmed by independent measurements. In the case of red 
giant stars recent work has cautioned about the predictions 
from scaling relations in first ascent red giant and red clump 


stars ( 
stars ( 

Miglio et al. 

20121 , expected masses o 

F metal-poor 
sroperties of 

Epstein et al. 

20141, and open cluster 

main-sequence binaries (Brogaard et al.||2012 Miglio et al. 


2012). Results from detailed modelling also suggest that the 


scaling relations underestimate the mean stellar density of 
red giant stars by ~5% ( Huber et al.|2013a K 

While acknowledging the shortcomings of our current 
understanding of the physics behind equations El and 
we note that the discrepancies arise in regimes of metal- 
licity and evolutionary state far from the main sequence for 
middle-aged stars near the mass of our Sun from where these 
relations are extrapolated. Our sample of KOIs presents an 
opportunity to compare determinations from individual fre¬ 
quencies with those of the grid-based method in a parameter 
space closer to the Sun, and thus to validate the results of 
the scaling relations for population studies in this region of 
the Hertzsprung-Russell diagram. The BASTA presented in 
Section [3.2| is flexible and can be adapted to any set of input 
data. We include as observational constraints the quantities 
Tefi, [Fe/H], (Au) and i^max of each target, which have been 


determined by Huber et al. (2013bI. To derive their stellar 
properties, we need to determine r'max and (Av) in stellar 
evolution models. 

The theoretical iZmax value can only be obtained using 
equation |12[ and we adopt for the reference solar value en¬ 
tering this equation i/max,© = 3090±30 /iHz given by|Huber| 


et al. (20111. The theoretical (Au) on the other hand can be 
obtained either from equation El or directly from the indi¬ 
vidual frequencies of oscillations calculated for each stellar 
model. To implement the latter approach we must take into 
account the way (Av) is estimated from the observed fre¬ 
quencies and treat models and data as similarly as possible 
(see White et al. I [MI] ). With this in mind, we adopt the 
following procedure: we determine the theoretical (Au) us¬ 
ing a Gaussian-weighted linear fit to the ^ = 0 modes of 
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Figure 14. Fractional differences in stellar properties deter¬ 
mined from the frequency combinations and the scaling rela¬ 
tions as a function of the large frequency separation, in the sense 
(combinations-scaling). 


each model as a function of radial order n, centred at the 
iZmax value obtained from equation and with a width of 
0.25 r^max. 

Figure shows the comparison between the stellar 
properties determined with the BASTA in section [^ using 
the frequency ratios and those obtained using the global 
asteroseismic parameters, where (Ai/) in the models has 
been computed using the aforementioned Gaussian fit to 
the theoretical frequencies. There are three outliers (same 
not shown in this hgure) corresponding 


4.1.2 


as in section 

to stars presenting bimodal PDFs. For the rest of the sample 
there is excellent agreement between both determinations, 
with median fractional differences and standard deviations 
of0.5±1.2% (density), 0.7±0.8% (radius), 1.8±2.1% (mass), 
and —10.6±10.8% (age). These results are reassuring consid¬ 
ering that median uncertainties in grid-based determinations 
are usually of the order of 2.2% in radius, 2.8% in density, 
5.4% in mass, and 25% in age ( [Chaplin et al.|[2014 |. Thus, 
stellar properties in the Tefi and [Fe/H] regime explored ob¬ 
tained from the global average asteroseismic parameters are 
compatible with those from individual frequency determina¬ 
tions and can be used when the latter are not available. 

Finally, we have computed a set of stellar properties 
for our KOI sample using the two scaling relations (equa- 
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tions 11 and 121 to determine the global asteroseismic pa¬ 
rameters, and the Huber et al. ( 2011| reference solar val¬ 
ues {{Avq) = 135.1 ± 0.1 fiRz, J/max,© = 3090 ± 30 /rHz). 
The median fractional differences and standard deviations 
with the results obtained from frequency combinations are 
—1.7±1.2% (density), 1.5±0.9% (radius), 2.6±2.2% (mass), 
and —10.4± 11.5% (age). The level of agreement is slightly 
worse than that obtained when the theoretical (Aiz) was 
computed from individual frequencies, but still well within 
the median uncertainties of the method quoted by |Chap-| 


lin et al. (2014l. However, compared to the results from 


the frequency ratios, the differences are of systematic na¬ 
ture with the scaling relations overestimating the density 
and underestimating the radius and mass. Departures from 
the {Au) scaling relation (equation |ll[ ) in stellar models can 
introduce biases in the results, and have already been re¬ 
ported by White et al.| poTI| and [Chaplin et ah (20141 in 
the main-sequence phase. A detailed exploration of this sys¬ 
tematic biases goes beyond the scope of this study and will 
be addressed in a subsequent publication. 


AGE CORRELATIONS WITH EXOPLANET 
PROPERTIES 


Using the stellar properties determined with the BASTA 
and shown in Table we can investigate the distribution 
of planetary properties such as orbital period, radius and 
multiplicity as a function of stellar age. Figure [^illustrates 
the age distribution of our sample as a function of planet 
radius and orbital period, showing that the majority of the 
high SNR asteroseismic host-stars are older than the Sun. 
Also, two main populations can be distinguished in age, at 
~3 Gyr and ~6 Gyr. These features arise due to biases in 
the detection of stellar pulsations, whose amplitudes are pro¬ 
portional to the stellar luminosity (e.g., |Kjeldsen fc Bedding 
1995). The net result is that detections are favoured in F- 


type (with ages about 2-3 Gyr) and old G-type (i.e., 6 Gyr) 
stars that have high enough luminosities (see e.g.. Fig. 1 in 
Silva Aguirre et al.||2011b l. Given these age biases for high 
signal-to-noise asteroseismic detections, it should not be in¬ 
ferred that our finding of an average age greater than solar 
applies to either the Kepler exoplanet host star sample in 
general, nor the larger set of all Kepler dwarf stars compris¬ 
ing the exoplanet search program. However, it is interesting 
to look for age correlations using this sample and explore 
the future prospects for larger samples of asteroseismically 
determined stellar properties of host stars. 

Figure shows no obvious trends of stellar age with 
planet radii or orbital periods. This is not surprising given 
that evolutionary effects which are expected to affect planet 
radii or orbital periods such as photoevaporation and orbital 


migration are expected to occur within <1 Gyr (Lin, Bo- 


denheimer fc Richardson|1996 Lopez fc Fortney|2013 1, and 

our sample does not contain such young stars. Nevertheless, 
our sample suggests that the period and radius distribution 
appears to be approximately constant for ages >2 Gyr, and 
hence is likely representative for Sun-like host stars observed 
by Kepler. We also do not detect any significant difference 
in the multiplicity fraction between young (~2-4 Gyr) and 
old (>5 Gyr) stars. 

Finally, Fig. shows the ratio of planetary periods 



Figure 15. Properties of single (red diamonds) and multiple 
(blue circles) systems as a function of stellar age. Top panel: plan¬ 
etary radius. Bottom panel: planetary period. The oldest multiple 
system correspond to the five planets detected around KOT3158. 
Vertical dashed line marks the solar. 


for all known multiple systems in our sample. Studies of the 
larger Kepler planet-candidate sample have shown that most 
planets are not on exact resonant orbits, but that their pe¬ 


riod ratios tend to be near first-order resonances (e.g., Fab- 
rycky et al.|2014 1 . Although our sample is small, the results 


seem to indicate that these conclusions hold independently 
of the age of the host stars. 


7 CONCLUSIONS 

We have made use of the full time-base of observations from 
the Kepler satellite to uniformly determine precise funda¬ 
mental stellar parameters, including ages, for a sample of 
exoplanet host-stars where high-quality asteroseismic data 
were available. We devised a Bayesian procedure flexible in 
its input and applied it to different grids of models to study 
systematics from input physics and extract statistically ro¬ 
bust properties for all stars. Our main results can be sum¬ 
marised as follows: 

The stellar properties determined from each grid of mod¬ 
els reflect the impact of the input physics in asteroseismic 
studies. The inclusion of microscopic diffusion clearly affects 




















































18 V. Silva Aguirre et al. 


14 

12 

10 
'XT 

^ D 

o 8 

(D 
(25 

< 6 

4 

2 

1.0 1.5 2.0 2.5 3.0 3.5 4.0 

PJP-. 

Figure 16. Period ratios for the multiple systems in our sam¬ 
ple. Dashed lines mark the 3:2, 5:3, 2:1, and 3:1 mean-motion 
resonances. 
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the ages derived while providing better hts to the spectro¬ 
scopically determined effective temperature and metallicity. 
Extending the convective core size by means of overshoot¬ 
ing results in slightly less massive stars than the standard 
Schwarzschild case once the convective core is well devel¬ 
oped and the target is close to the TAMS. Interestingly, we 
show that including this effect favours a particular solution 
when bimodal distributions in mass are present. 

Our results show that statistical uncertainties are below 
2% in radius and density, ~3% in mass, ~4.5% in distance 
and ~14% in age. This level of precision is perfectly suited 
for supporting characterisation of exoplanet systems, includ¬ 
ing determination of eccentricities by comparison of stellar 
densities from transit measurements with asteroseismic ones 


(e.g., Fogtmann-Schulz et al.|[2014 Van Eylen et al.]|2014 |. 
Furthermore, our catalogue of stellar ages can be used for 
studying time evolution of planetary systems, tidal interac¬ 
tions and circularisation timescales. 

We have tested the systematic uncertainties introduced 
by changing the solar abundances and the efficiency of con¬ 
vection, and found them to be smaller than the statistical 
errors. We explored the impact on the stellar properties of 
reproducing different sets of asteroseismic observables (indi¬ 
vidual frequencies or frequency combinations), using differ¬ 
ent evolutionary and pulsation codes, as well as fitting algo¬ 
rithms, and found that when the input physics is kept the 
same the systematic differences are of the order of 1% in den¬ 
sity and radius, and 2% and 9% in mass and age. A source of 
uncertainty comparable to the statistical ones comes from 
the relation between metallicity and initial helium (Vni), 
and we expect the Gaia results to produce determinations 
of luminosity and provide an independent constraint on the 
initial helium abundances. Fitting algorithms that do not 
restrict this parameter commonly predict too low values of 
the Yini, even below standard big bang nucleosynthesis. In 
Kepler-21, where a parallax measurement is available, the 
obtained luminosity of the sub-SBBN model is 1.5 cr away 
from the parallax inferred value. 

When the SNR of the data is too low and individual fre¬ 


quencies of oscillation are not possible to extract, astero¬ 
seismic determination of stellar properties is made using 
the global asteroseismic parameters. We have applied the 
BASTA using these observables and the spectroscopic con¬ 
straints as input and found an excellent agreement with 
properties determined from frequency combinations. These 
results validate the use of the global asteroseismic observ¬ 
ables in the regime explored in this paper when the theoreti¬ 
cal value of (Al'} is computed from individual frequencies of 
oscillations. If determined from the scaling relation instead, 
the results should be treated with caution as they seem to 
overestimate the density and underestimate the mass and 
radius. 

Our results further demonstrate the positive synergy be¬ 
tween asteroseismology and exoplanets studies, with the 
added ingredient of precise age determinations to a level 
better than 15%. We have investigated age correlations with 
planetary period and radii in our sample, finding no notice¬ 
able trend in the results. Their distribution as a function of 
age appears constant and is likely representative of the un¬ 
derlying population of Kepler Sun-like host stars. In terms 
of resonances, our findings suggest that the period ratios 
in multiple systems tend to be near first order resonances 
regardless of the age of the system. Although subject to se¬ 
lection biases, our sample shows that the majority of planet 
host stars with asteroseismic observations of high SNR are 
older than the Sun. 

The results presented here will provide the basis of up¬ 
coming catalogues of asteroseismic properties of exoplanet 
host-stars, where over 100 KOIs with (Au) and r'max mea¬ 
surements will be characterised (Huber et al. 2015, in prepa¬ 
ration). This larger sample will extend the investigation we 
have initiated in this paper on the dependance of plane¬ 
tary properties with the age of the host star. Similarly, 
the BASTA combined with grids of evolutionary models 
will provide stellar properties for targets observed by the 
TESS and PLATO missions in evolutionary stages beyond 
the main-sequence phase. 
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Table 3. Recommended set of stellar properties and statistical uncertainties of exoplanet candidate host stars as determined with BASTA. Solar luminosity used Lq = 3.846x 10^^ (erg/s). 
Additional systematic uncertainties can be accounted for (see sections |4.2j[X^ and |4.4] : from input physics of 0.8% (density), 0.7% (radius), 2.3% (mass), 9.6% (age); from choice 
of observables 0.3% (density and radius), 1% (mass), and 7% (age); from fitting algorithms and codes of 1% (density and radius), 2% (mass), and 9% (age); from the initial helium 
abundance of 1.7% (density), 1.6% (radius), 3.6% (mass), and 16.8% (age). Reference column gives literature sources of confirmed or validated exoplanets, ^see also |Fressin et aT~| | |201l] . 
'^Ephemeris match indicates contamination, see CFOP website (https://cfop.ipac.caltech.edu/home/) 
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APPENDIX A: EVOLUTIONARY CODES AND 
FITTING ALGORITHMS 

In the following sections we describe the fitting procedures used 
to determine stellar properties for our sample, including the evo¬ 
lutionary and pulsation codes used. 

Al ASTFIT 

The ASTEC Fitting method (ASTFIT) uses the Aarhus stellar 
evolution (ASTEC, |Christensen-Dalsgaard|2008b[ | and pulsation 
(ADIPLS, |Christensen-Dalsgaard|2008a[ | codes. It is based on fit¬ 
ting individual model frequencies, applying a correction for the 
surface errors, to the observed frequencies (see also |Gilliland et ah| 
|2013[ |. The fit is carried out in grid of models varying the mass 
and composition. The chemical composition is based on a he¬ 
lium enrichment relation AYfAZ = 1.4, set up in terms of ini¬ 
tial hydrogen and heavy-element abundances (Aq, Zq), with three 
neighbouring values of Zq for each value of Aq. The grid extends 
in [Fe/H] between —0.577 and 0.415 in steps of around 0.05 and 
in mass between 0.7 and 1.69 Mq, in steps of 0.01 Mq. The 
ASTFIT results presented here were computed without diffusion 
and settling, and convective-core overshoot was not considered. 
The mixing-length parameter was set to cimlt = roughly 
corresponding to solar calibration. The model physics otherwise 
essentially followed the BASTA specifications. 

The fit of a given model to the data is defined in terms of 



where and are the observed and model frequencies 

and (Ji is the error in the observed frequencies. Here the model 
frequencies typically include a surface correction (see section [Al.l| 
below). It is assumed that the degree and order of the observed 
frequencies have already been determined, and that they do not 
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include additional mixed modes. Hence, when the models con¬ 
tain mixed modes, they are excluded (renumbering the remaining 
modes) retaining only the most acoustically dominated modes (as 
determined from their inertia) in each interval in large frequency 
separation Au. The fit also includes 


^jobs) _ ^jmod) y ^ 

ATes) y <^([Fe/H]) 

with the observed and model effective temperature Tgff and 
[Fe/H]. The complete fit is then characterised by = x5 + Xspec? 
defining the relative weight between the oscillation and the spec¬ 
troscopic data. 

For each evolution sequence frequencies are calculated for 
at least every second timestep. Amongst these frequency sets 
the model A4' . with the smallest y?, is determined. Based 

min 

on homology scaling it is assumed that the frequencies in the 
vicinity of obtained as where r = 

R being the surface radius of the model, and 
the best-fitting model is determined by minimising 



as a function of r. The resulting value Tmin of ^ defines an esti¬ 
mate -Rmin of the radius of the best-fitting model along the given 
sequence. 

The other properties of this best-fitting model, for the given 
set of model parameters {Vk} with corresponding 
are determined by linear interpolation in R to -Rmin- The corre¬ 
sponding minimum . {Vk) is then determined as y^ . {Vk) = 
X^,inin(^*) + Xapec(Vi;jjin|^ The final best-fitting model is found 
as the parameter set corresponding to the smallest X^niin^^^) 
(or Xmin(^^)) over a suitable selection of models in the grid. 

It is very informative to plot y^ min(^^) against 

the various model properties, such as mean density (p), R or 
age. Representative results are given in the HAT-P-7 paper 
< |Christensen-Dalsgaard et al.|2010| |. The best-fitting frequencies, 
e.g., for comparison with the observations in an echelle diagram, 
are obtained by applying the appropriate scaling rmin to the fre¬ 
quencies of the model in the minimising sequence. 

Although probably not completely justified in a statistical 
sense, estimates of stellar properties, in particular radius, mass 
and age, are determined as averages and standard deviations 
of the properties of the models AiminiRk) over the parameters 
{'Pk}'> with the weights exp(—y^.j^(7^;c)/2), corresponding to the 
likelihood. These are the final ASTFIT results compared with the 
BASTA results in figures and \TT\ and presented in Tabl( jAl] 




A 1.1 Near-surface correction 

It is well established that errors in the modelling of the struc¬ 
ture of the near-surface layers of the star and their effects on the 
oscillations cause systematic differences between the model and 
observed frequencies, even when the model otherwise matches the 
structure of the star. This is particularly evident in the solar case, 
where the availability of frequencies over a broad range of degrees 
allow separation between the effects of the bulk of the solar in¬ 
terior and the surface layers. These errors must be taken into 
account when, as done in ASTFIT, individual model frequencies 
are fitted to the observations. The correction unavoidably requires 
prior and perhaps poorly justified assumptions, typically based on 
the behaviour determined in the solar case. 

^ and hence neglecting the generally negligible change in, e.g., 
Teff between timesteps. 



2000 3000 4000 5000 

V (//Hz) 


Figure Al. The solar surface correction determined from 

7^2 in the fit given in Eq. ( |A4[ | to solar data over a broad range 
of degrees. 


Here we specifically use a representation of the surface 
correction that has a functional form matching the solar be- 
haviour (see also |Aerts, Christensen-Dalsgaard &: Kurtz| 201Q[ 


|Christensen-Dalsgaard|2012[ |. This is obtained from the differen- 

tial form of the Duvall asymptotic expression for the frequencies 
( [Christensen-Dalsgaard, Thompson &: Gough|1989| l, according to 
which the frequency differences Su between the Sun and a model 
satisfy 

Snl - — — R^iynllN) -h ^ (A4) 

^nl 

where L = -\Jl(l -h 1) and Sni is a scale factor which can be 
calculated from the model, and which may be chosen to be close 
to 1 for low-degree modes. The term in Ri reflects the difference 
in sound speed between the Sun and the model throughout the 
star, whereas 7^2 contains the contribution from the near-surface 
region. Since on general grounds the surface effects are very small 
at low frequency, the arbitrary additive constant in the definition 
of 7 ^2 m ay be chosen such that 7^2 is zero at low frequency. In 
is shown the frequency shift caused by the surface 


Al 


effects in the Sun, determined in this manner from 1/7^2 (^)- We 
now make the assumption that the correction in other stars has a 
similar functional form, in terms of frequency measured in units 
the acoustic cut-off frequency which provides a physically 
motivated frequency scale for the stellar atmosphere. Thus for 
each set of model frequencies we obtain the surface correction by 
determining a scale factor r and an amplitude d from a least- 
squares fit of 


I^nl = + agQiunilus.c) 


(A5) 


to the observed frequencies, where the function Qq is determined 
from illustrated in Fig. using a solar acoustic cut-off 

frequency of z^ac,© = 5200 /zHz. The model frequencies in equa¬ 
tion are obtained from this expression, but without including 
the scale factor. 


A2 YMCM 

The YALE Monte-Carlo Method used the YREC ( jPemarquej 
jet al.|2008| l code to model the stars. While the input physics was 
the same as in the other techniques, the method of constructing 
the models was very different. 

The starting point of the modelling effort for each star was 
the average large separation that was used in the grid-based mod- 
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Table Al. Stellar properties determined with ASTFIT. Solar luminosity used Lq = 3.846 x 10^^ (erg/s). The first five objects are 
shown here for guidance on the format, and the full table is available online. 


KOI 

KIC 

Mass 

(Mq) 

Radius 

(-R©) 

Density 

(g/cm^) 

log g 
(dex) 

Luminosity 

(R©) 

Age 

(Gyr) 

2 

10666592 

1.538ib 0.030 

2 .000± 0.012 

0.271± 0.002 

4.023± 0.004 

5.S11± 0.213 

1.88± 0.14 

5 

8554498 

1.173± 0.020 

1.785± 0.010 

0.291± 0.001 

4.004± 0.003 

3.373± 0.133 

6 .12± 0.28 

7 

11853905 

1.118± 0.031 

1.553± 0.013 

0.420± 0.002 

4.104± 0.005 

2.488± 0.117 

7.27± 0.69 

41 

6521045 

1 .121± 0.012 

1.517± 0.006 

0.452± 0.001 

4.125± 0.002 

2.509± 0.088 

6.83± 0.28 

42 

8866102 

1.218ib 0.036 

1.354± 0.013 

0.691± 0.004 

4.260± 0.005 

2.654ib 0.111 

2.78± 0.47 


elling, the mass obtained from grid-based modelling and the ef¬ 
fective temperature and metallicity obtained from spectroscopy. 
Since each of these quantities is associated with an uncertainty, 
we created many more realisations of these parameters to obtain a 
larger set of (M, Au, Tgff, [Fe/H]). The Av in each set, along with 
M was used to calculate R using the usual scaling relationship. 
Since the scaling relationship is not exact, we used the empirical 
correction of |White et al.| | |2011[ | to get a better estimate. 

For each combination (M, Teff, [Fe/H]) we used YREC 
in an iterative mode to obtain a model of the given mass M and 
[Fe/H] than had the required R and Tgff. Since the mixing length 
parameter was set to be the solar calibrated value (1.6756 for the 
given physics and formulation of mixing-length theory in YREC), 
the only free parameter we had was the initial helium abundance 
Yq that was varied until the required model was obtained. Some 
combinations of parameters resulted in models that need Yq to 
be lower than the primordial helium abundance; all such models 
were rejected. 

Theoretical freque ncies of oscillat i on are computed using the 
code first described by |Antia &: Basu] ( |l994^ . The for fho dif¬ 
ferences between the observed frequencies and those of the mod¬ 
els, were calculated after correcting for the surface term in 
a manner outlined below. Similarly, we calculated Xrat f®*^ 
differences in the frequency ratios roio and ro 2 (including cor¬ 
relations), Xt ff f*^^ difference in effective temperature, and 
^fpe/H] difference in metallicity. The final goodness-of-fit 

was determined by 

Xtot = Xi + Xrat + + Xfpe/H] ' 

The best-fit model was the one with the lowest value of Xtof 

We used the scaled solar surface term to correct the fre¬ 
quency differences between the observed and model frequencies. 
We first constructed a solar model with the same physics as the 
rest of the models and determined its frequency differences with 
respect to observed solar frequencies obtained by the Birming¬ 
ham Solar Oscillations Network (BiSON) and listed in Table 1 of 
[Chaplin et al.| ( [2007| l. The solar surface term was obtained by fit¬ 
ting these frequency differences to the differential form of the Du¬ 
vall Law ( jChristensen-Dalsgaard, Thompson & Gough|1989| |. We 
denote this as the © relation. Both © and Su^i q 

are then scaled to the mass and radius of the stellar model under 
consideration using the homology scaling 


(Ai^(mod)) 


(A7) 


where the angular brackets denote the average. The resulting 
'^^ni © is then used to correct the stellar model for the 

surface term. The corrected frequencies are denoted as 


fcorr) (mod) 

nl nl 


gr&Unl,@, 


(AS) 


with r5u„i q evaluated at ru„i 0 =i/„;(obs). The factor /3 is se- 
lected to minimize 




(obs) (corr) ' 

j ' 

nl nl 




(A9) 


are then used to calculate xj- The final set of stellar prop- 
erties determined with the YMCM are available in Table [A2] 


A3 AMP 

For the results presented in s ection |4.3| we used the Asteroseis- 
mic Modeling Portal (AMP, [Metcalfe, Creevey &: Christensen -j 
[Dalsgaard|2009] [Woitaszek, Metcalfe &; Shorrock|2009[| in a sim¬ 

ilar configuration as that described in [Metcalfe et al.[ ( [2014[ |, but 
without including diffusion and settling of helium. In summary, 
AMP uses a parallel genetic algorithm (GA, [Metcalfe & Charbon-[ 
[neau[[2003| l to optimize the match between stellar model output 
and the available set of observational constraints. The evolution 
models are produced with ASTEC, and the oscillation frequencies 
are calculated with ADIPLS. The five adjustable model parame¬ 
ters include the mass (M), age (t), composition {Z and Yi), and 
mixing-length (a). The oscillation frequencies and other prop¬ 
erties of each model are compared to four sets of observational 
constraints, including: [ 1 ] the individual frequencies corrected for 
surface effects following the empirical prescription of [Kjeldsen,[ 
[Bedding Christensen-Dalsgaard[ ( [2008[ l, [2] the frequency ratios 
r 02 defined by Eq.[^ [3] the frequency ratios roio defined by Eq.[3 
and [4] the available spectroscopic constraints. A normalized ^ 
is calculated for each set of constraints, and the GA attempts to 
minimize the mean of the four x^ values. This allows the various 
asteroseismic quality metrics to be traded off against each other, 
while ensuring that the numerous frequencies and ratios do not 
overwhelm the relatively few spectroscopic constraints. The in¬ 
terested reader can find the optimal model parameters produced 
by AMP for each star in Table [A^ 
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Table A2. Stellar properties determined with YMCM. Solar luminosity used Lq = 3.846 x 10^^ (erg/s). The first five objects are shown 
here for guidance on the format, and the full table is available online. 


KOI 

KIC 

Mass 

(Mq) 

Radius 

(«0) 

Density 

(g/cm^) 

log 9 

(dex) 

Luminosity 

(L@) 

Age 

(Gyr) 

■^ini 

Zini 

2 

10666592 

1.546± 0.035 

2 .001± 0.016 

0.272± 0.002 

4.024± 0.004 

5.744± 0.175 

1.63± 0.10 

0.278±0.017 

0.0282±0.0038 

7 

11853905 

1.146± 0.030 

1.567± 0.014 

0.419± 0.005 

4.107± 0.005 

2.533± 0.100 

6.93± 0.55 

0.257±0.013 

0.0232±0.0026 

41 

6521045 

1.138± 0.009 

1.525± 0.004 

0.452± 0.001 

4.127± 0.001 

2.564± 0.048 

6.22± 0.14 

0.267±0.006 

0.0219±0.0020 

42 

8866102 

1.238± 0.041 

1.359± 0.015 

0.694± 0.003 

4.264± 0.005 

2.671± 0.083 

2.32± 0.31 

0.268±0.016 

0.0187±0.0021 

69 

3544595 

0.919± 0.020 

0.922± 0.007 

1.650± 0.002 

4.471± 0.003 

0.771± 0.028 

7.04± 0.68 

0.252±0.015 

0.0133±0.0015 


Table A3. Stellar properties determined with AMP. Solar luminosity used Lq = 3.846 x 10^^ (erg/s). The first five objects are shown 
here for guidance on the format, and the full table is available online. 


KOI 

KIC 

Mass 

(Mq) 

Radius 

(«©) 

Density 

(g/cm^) 

log g 
(dex) 

Luminosity 

(L@) 

Age 

(Gyr) 

yini 

Zini 

Q^MLT 

2 

10666592 

1.300± 0.055 

1.879± 0.018 

0.276± 0.001 

4.004± 0.003 

5.096± 0.291 

2.79± 0.52 

0.317±0.011 

0.0239±0.0043 

1.72±0.13 

5 

8554498 

1.310± 0.038 

1.825± 0.018 

0.303± 0.001 

4.032± 0.005 

3.444± 0.190 

4.56± 0.60 

0.243±0.011 

0.0255±0.0032 

1.44±0.09 

7 

11853905 

1.140± 0.017 

1.553± 0.016 

0.428± 0.008 

4.112± 0.007 

2.435± 0.185 

6.13± 0.30 

0.268±0.012 

0.0217±0.0027 

1.62±0.12 

41 

6521045 

1.090± 0.018 

1.501± 0.013 

0.454± 0.008 

4.123± 0.005 

2.337± 0.129 

7.16± 0.42 

0.265±0.010 

0.0184±0.0012 

1.78±0.07 

42 

8866102 

1.260± 0.032 

1.363± 0.014 

0.701± 0.004 

4.269± 0.003 

2.626± 0.137 

2.01± 0.47 

0.257±0.012 

0.0184±0.0015 

1.82±0.09 



